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l.O INTRODUCTION 


The work done by ACM under Contract NAS9-15171 covered 
a wide variety of tasks. The categories of work includ('d: 

• Development of ».ew mathematical and al gor 1 1 hmi cal 
techniques for solution of jiroblems In satellite 
dynamics . 

• Development of new solutions to satellite orbit 
mot ion . 

• Dynamical analysis of Shuttle on-orblt operations. 

• Development and analysis of computer software routines 
for use in Shuttle mission planning. 

• Formulation of mathematical models of atmo.spherlc 
density . 

Due to the variety of the tasks and the several modifi- 
cations to the contract, a separate report has been written 
for each task. These reports are listed in the references. 

The purpose of this Final Report is to describe the objec- 
tives and scope of each task, and to give a summary of the 
significant results of each task. Therefore, this report 
will serve as an overview and summary of the work done under 
Contract NAS9-15171. In case more detailed information is 
required, please refer to the individual reports. Copies of 
each report are available from the Technical Monitors or the 
JSC Technical Library. 
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2.0 ANALYTICAL SOLUTION OF ARTIFICIAL 
SATKI.LITE MOTION WITH DRAG 

Orbit computation methods can usually be classified 

as : 

Num erical - The calculations are carried out in a step- 
by-step manner. High precision is i>osslble but 
computer runtime can be excessive. 

An al yt ical - The calculations are carrlt*d out in one 

step, regardless of the prediction interval. There- 
fore, these methods liave extr(.‘mely fast computation 
t Imes . 

The fast execution tlmc»s of analytical satellites theories 
make them very attractive as mission analysis and planning 
tools. But difficult problems existt»d in the analytical 
met hods : 

1) The solutions were expressed by extremely complex 
formulas which recpiired much more computer storage than the 
numerical methods. 

2) The accuracy of the analytical methods was restricted 
by the simplicity of the model or the solution technicjue. 

3) Unlike numerical methods, the analytical theories 
were generally difficult to extend. For instance, numerical 
methods may use a goopotential model of arbitrary order or 
degree, where as the analytical methods wore limited to a fixed 
order and degree. The numerical method.s can easily include 
new or more accurate density models but this was impossible in 
analytical methods. 

4) Analytical theories suffered because the different 
perturbations had never been unified under one non-singular 
formulation . 

These problems were the motivation behind developing a 
new satellite theory. The resulting theory provides a con- 
cise and efficient algorithm which has been implemented into 
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a working, modular computer software tool. The development 
of the theory and corresponding software are described in the 
following sections. 


2^1 Perturbat icms due to AlnK)y)herlc Draj[ 

The objective of this task vas to arrive at a fully ana- 
lytical theory for the motion of an artificial satellite per- 
turbed by the Jj term of the zonal geopotential expansion 

and by a drag force which is tangential to the orbit and pro- 
portional to the square of the velocity magnitude. The task 
was executed in two stages: 

1) Development of the drag differential equations in a 
regular set of two body elements and in a form suitable for 
solution by quadrature (reference 1). 

2) Development of a realistic density rm»del suitable 

for application in the theory developed in stage 1 (reference 7). 

The second stage is discussed in more detail in Section 4.2 
of this report and so this section will restrict its discussion 
to the first stage. 

Thr drag theory was built on top of the previously devel- 
oped theory in the PS^ elements. Since the density ef- 

fects only the magnitude, not the direction, of the drag force, 
the transformation of the drag forces may be carried through 
without specifying the density function. A canonical treatment 
of this transformation of the drag effects results in a tract- 
able set of differential equations. To allow the analytical 
solution of these equations, an expansion is made in powers 
of the eccentricity. To avoid errors, these expansions were 
created in an automatic manner by a minicomputer program spe- 
cially designed for this task. The result is a theory of rela- 
tively concise form, programmed in FORTRAN and executable in 
the interactive mode on the UNIVAC 1110. 

Numerical studies were performed comparing the analytical 
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Holution to a precision numerical solution. In the first step, 
both analytical and numerical solutions assumed the density 
to be a constant. The results showed the analytical constant 
density solution to be accurate to within of the numerical 

solution for Shuttle type orbits. In the second step, the 
analytical theory with the realistic density model was compared 
to num<*rlcal solutions obtained by using <*xtremely accurate 
and complex drag models. Kvon In those rigorous tests, the 
analytical solution showed errors of no more than 2^ 

2^2 Pertu rba t ions due to an Arbitrary Number of Zon al T^rm^ 

It was the purpose of this task to derive a "complete" 
first order solution to the motion of a satellite about a ro- 
tationally symmetric earth. No restriction was made on the 
order of the geopotential expansion, and the solution was 
developed entirely in non sit.gular elements. Recursive ex- 
pressions are used so that a concise formulation is maintained. 

The method of solution was tc apply the Von Zeipel al- 
gorithm to a recursive formulation of the zonal geopotential. 
This is done by defining the zonal potential in the canonical 
and regular PS^ elements. Considerable simpl icat ions can 
be realized with this choice of elements. By including the 
true anomaly as an canonical element, the hamiltonlan for a 
finite number of zonal harmonics becomes a finite fourier ex- 
pansion in the canonical elements. Also, the second order 
integration of the mean motion is not necessary because the 
new formulation uses a mean motion related to the total energy. 
The detailed development of this theory is given in reference 2. 

The result is a simple formulation which has been implemented 
in FORTRAN and executable on the UNI VAC 1110. 

Numerical studies were conducted to determine the accuracy 
of the solution. In the studies, the analytical solution was 
compared to a precision numerical solution. Both numerical 
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Hnd analytical Holutions modoled the IRth order geopotential 
expansion. For Shuttle-type orbits, the differences betw<*en 
analytical and numerical solutions were on the order of a fewr 
meters after 500 revolutions. 

2^3 An Analytical Satellite Orbit Pred ic tion Program (ASOP)^ 

The analytical formulations that are discussed In Sections 
2.1, 2.2 and 4.2 of this report as well as certain formulations 
that were developed under a previous contract, are being con- 
verted Into an operational computer program. This program rep- 
resents a new state-of-the-art for orbit prediction computer 
software. ASOP is able to compute near-earth orbits to within 
an accuracy of a few meters. Recursive equations are used 
instead of complicated formulas so that the execution time of 
ASOP Is on the order of a few milliseconds. At present those 
calculations only Include the perturbations due to Jj . How- 
ever, the perturbations due to atmospheric drag and all the 
Earth’s geopotential terms are currc^ntly being incorporated 
into the existing ASOP program. 

The current version of ASOP le now available to the 
general user of the UNIVAC 1110-EXEC 8 system at NASA/JSC 
Houston in two forms: 

(1) A stand-alone program that can be used Interactively 
to obtain Immediate results for a specific problem. 

(2) A user-subroutine package that can be Incorporated 
into other software systems. 

Both versions were designed to be small and execute quickly. 

ASOP has been extensively documented (reference 3). Any 
additions to the program will be documented as revisions to 
reference 3. 
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3.0 AN ANALYTICAL STATE TRANSITION MATRIX FOR ORBITS 
PERTURBED BY AN OBIJn'E SPHEROID 

Oftt<n In orbit dett>rmi nat Ion or naviKutlon algorlthma 
ont* muMl prodict how small dovtationB from some nominal orbit 
will in lime causo the actual orbital path to d«>vlate from the 
nominal path. If the initial deviation is small, then a 
linear approximation may be used to detc*rmine the deviation 
at any ^iven time. The linear approximation, a truneat<*d 
Taylor series t>xpanKion, requires the first partial deriva- 
tives of the cotji’dinates defininf? the position and velocity 
at a given time, with respect to the coordinates defining the 
position and velocity at the initial time. This matrix of 
partial derivatives is commonly called the state transition 
matrix - ♦ 

may be determined by numerical Integration. But when 
one considers the fact that the transition matrix is usually 
used in some iterative manner, the computational costs re- 
quired by the numerical Integration of the matrix differential 
equation becomes prohibitive. 

The alternative, of course, is to determine the transition 
matrix by some analytical technique. Several analytical meth- 
ods have determined expressions for the transition matrix under 
the assumption that the satellite moves along a two-body orbit, 
all perturbations being neglected. But this assumption may 
result in non-negl igible errors if the perturbations are large. 

For artificial satellites orbiting near the earth, the 
oblateness potential contributes a strong perturbation 

which may not be neglected for accurate satellite orbit pre- 
diction. Significant advantages are offered by a satel- 

lite theory developed from a new set of canonical €‘lements 
proposed by Schelfele. The elements are in an extended phase 
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Kpace in that eight (Instead of six) variables describe the 
state. They also have as their independ(>nt variable, a quan- 
tity related to the true anomaly Instead of time. The set is 
similar to the Poincare elements and are therefore named the 
Poincare-Similar (PS) elerm^nts. The PS pli*ments are non- 
singular for bounded orbits. 

Because of the simplicity and accuracy of the canonical 
analytical solution bused on PS elements, it becomes c log- 
ical choice on which to base a new analytical state transition 
matrix. Therefore the intention of this task was to derive a 
completely analytical singularity free form of the state tran- 
sition matrix for orbits perturbed by an oblate spheroid. 

The result is a concise, accurate and algorithmically 
simple state transition matrix and its Inverse which includes 
the short period and secular effects of the xonal harmonic. 

The algorithm has been Implemented in a working program exec- 
utable in interactive mode on the UNIVAC 1110. From numerical 
studies it was found that the improvement over the conic matrix 
ranges from 2 to 4 digits better accuracy. Complete docu- 
mentation on the work of this task is given in reference 4. 
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•i.O DKVKI.OPMKNT AND APPLICATION OF 
ATMOSPHKRIC DKNSITY MODKI^ 

For uccurnte orbit oomputut lonH fur HutollitoH porturbod 
by atmosphorlc drag, h rathor complc*x dopHlty model must bo 
<*mployed. The Jacchia density m»)del includes several variations 
of the atmosphere character which are related directly or indi- 
rectly to the density. In the most recent Jacchia model, the 
exospheric tc*mperature .is modeled by variations due tot 

1. Averaged and dally variations in the solar flux (*^j0 7^ 

2. Averag(>d and 3 hourly variations in the geomagnetic 

index (a ) 

P 

3. Variations of the angle betwret'n satellite and atmos- 
pheric bulge (diurnal effect) 

The exospheric temperature is then relat<*d to the density by 
the solution of the diffuse »*qullibrium equation for the dif- 
ferent constituents of the atmosphere and the desired altitude. 
Other variations are then computed directly as changes in the 
density; the> include; 

1. Changes due to st*mi annual effect. 

2. Changes due to seasonal-latitudinal effect. 

The causes for these variations are not exactly known but may 
be modeled by fitting emperical data. 

Although the mf*del is quite accurate, it requires a great 
deal of computation time and storage. The Inefficiency of the 
model can not be reccnclled by simply neglecting some of the 
effects and other models such as the 1962 standard are grossly 
inaccurate for altitudes above 90km . In addition, the com- 
plexity of the Jacchia model renders it useless for analytical 
satellite theory applications. 

The goal of this task was therefore to develop new density 
models which are efficient for applications in numerical methods 
or which may be applied in satellite theories. The manner in 
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which thc*HC KonlH were met are descrlhed in detail in refer- 
encoH 5, 6 and 7. A brief overview in Klven in the followinK 
Hee t lt)n« . 


A.l_ Appll cat ion to Numerical Orbit Computation Methods 

In an effort to avert the computational costs, two density 
models have been Independently develt)p«»d to slmulato the Jacchla 
iTKidel without the costs of an Inefficient alKorithm. The first 
model, referred to as nubb-Muellur (B-M), has bt*en develop<*d 
In Mueller (reference 5) and numerical studies showed it to 
be an efflcti*nt yot accurate model. The B-M model assumes the 
density can be modeled by rather simple expresKions but the 
parameters of the mctdel are unknown. These c^xpressions can 
be easily Inverted l.e. Riven the di'nslty at a few points the 
parameters can be determined. Thus the B-M model can be "cal- 
ibrated" to fit the Jacchla model, by evaluatinR the Jacchla 
mod€»l at a few specific points. Tlie model parameterr found 
by calibration are Implicit functions of the slowly varying 
characters of the atmosphere and can be considered constants 
over a limited pt*rlod of about a month. 

The second mode*l , n:uned Jacchia-Mnebc*rry (J-L), is of 
a similar nature as the B-M model. It also assumes rather 
simple expressions for the density whose paramc'ters have been 
determined by a least square fit with the Jacchla model. Un- 
like the B-M model, the J-L model’s parameters are indeed con- 
stants and thus no maintenance is required to update the para- 
meters. But the J-t model does have an overhead in that a 
considerable number of parameters are required for the com- 
plete model. In any case, both models show large computer cost 
savings over the complex Jacchla model. 

In this task the algorithm for the B-M model was documented 
in reference 6. Also in that reference, the complete develop- 
ment of the J-L model was given. Sevt*ral versions of the J-L 
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niodol were d€^volop»*d iind Impl «*d tn n r<»mpiiti»r parka^o on 
the UNIVAC 1110. A usor'n 1« p:iv«-n In rc^fcT»*nce 6, Inrlud- 

Ing Hoftware dorumont at Ion . In addition, niimorlcal Hludlefc Hhow»*d 
the J-Ii model to be very accurate for altltudf'M bolow 500km 
Hun t ImoH ahowed the J-I. model to be 3 to 9 timoK fnater 
than the Jacchla model depending on which veralon wna UK<>d. 

•1.2 Application to Analytical Satellite Thcorlo^ 

In the puHt, the computation of the drag p<*rturbat ions 
by analytical methods has not been feasible because of their 
inability to model the "real world" density. 

In developing a density model for the analytical theory, 
one Is sev<*rely restricted by the fact that the model must be 
In the form of a fourier series In the true longitude. As Is 
the case In most analytical theories, the perturbation must 
be written in a ff)urier series to facilitate the solution of 
the differential e(]uat ions of motion. Usually, the soluti{)n 
is obtained by quadrature. 

Several density models have been dt‘ve loped to predict 
accurately the density at any point in space and time. Exam- 
ples are the Jacchia model and the USSR model. But both models 
are extremely complicated and too unwieldly for analytical 
satellite theories. In the analytical theory of Brouwer and 
Hori, the density was assum<^d to be an exponential function 
of the position radius of the satellite. The exponential must 
then be expanded in a Poisson series so the quadrature can bo 
performed. This model has several difficulties. It first has 
a problem with convergence, which Brouwer points out. Secondly, 
it is simply a poor model for describing the dynamic atmosphere. 
The density is extremely effected by such factors as the level 
of solar activity and whether it is summer or winter, day or 
night. Thus the model in Brouwer. Hori theory is simply in- 
adequate . 
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RoCv*ntly an extrfmoly simple density model (referred to 
here as B-M) has l)een developed to match the Jacchia model to 
a high dejfroe of accuracy (reference 5). The variations in 
the density due to changes in the height and changes in the 
relationship of the vehicle and sun position (diurnal effect) 
arc includtfd explicitly. I.ong period variations such as the 
changes in the average solar activity and semi-annual vari- 
ations are included implicitly in the coefficients of the model. 
The value of the coefficients are determliu'd by a procedure 
called "calibration". The simple formulation allows the model 
to be inverted, i.e. glvc'n the density at different points in 
space (as determined from Jacchia) one can compute the coeffi- 
cients of the B-M model. Since the coefficients are implicit 
functions of long period effects, they can be considered con- 
stants over a limited period of about a month. 

Fven this extremely simple model cannot be appli(‘d in the 
analytical drag theory because it cannot be written in the form 
of a fourier series. However the technique of the B-M model 
does give important insight and direction to follow. 

In all the models discussed, the represontat '.ons of the 
density are considered to bo gloi)al . In other words, given 
any position In the atmosphere one can determine its density. 

The approach taken in this task was to develop a modc'l which 
expresses the dens'* *.y along a particular orbit. The coeffi- 
cients in the model are calibrated with the Jacchia model in 
u manner similar to the B-M model. But in this case, the co- 
efficients in the model are not only implicit functions of the 
long period variations in the atmosphere but also the orbital 
elements which describe the orbit. The result is a density 
model which can be written in a fourier series and easily im- 
plemented into the drag theory. Since the orbital elements 
are perturbed by and drag, the coefficients in the model 

must be updated periodically or corrected in some manner to 

If 

F 
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ref lect the changes. Since the perturbations are small, the 
updating would bo infrequent. DiMails of the devolojiment can 
be found in reference 7. 

The new density model has several distinct ndvaulages. 

It has a rather subtle yet c»xtremely important advantage in 
that the model is a power expansion in Cj I l.e., the basic 

expansion used in Seheifele's development of the dl f f er<*nt 1 al 
(>(|uat ions. 

Another advantage of the new density model is its ability 
to simulate any density model. In this task the Jacchia model 
was chosen for the simulations, but any other model could have 
been used. If more accurate density models are developed, they 
may also be chosc'n for simulation and thus the analytical drag 
theory can reflect the accuracy of the newly developed models. 

Finally, it is concluded that the analytical satellite 
theory need not be limit(‘d by a simplified drag model. With 
the approaches developed in this report, it is feasible to 
make the accuracy of the analytical theory competitive with 
preoieion numerical methods, while retaining a concise formu- 
lation and quickness of execution. 




-19- 


5.0 OHHITAL McmON OK THE SOI-AR POWER SATEM.ITE 

It has been proposed to put a series of large satellites 
Into geosynchronous orbit for the purpose of collecting solar 
energy and redirecting It toward the <>arth via microwave radi- 
ation. Preliminary studies are b(*ing carried out at JSC on 
the feasibility of these Solar Power Satellites (SPS). 

Tlie large area of the collecting surface (approximately 
Idd square kilometers) means that so''ar radiation pressure 
will cause significant perturbations on the SPS orbit. In 
fact solar pressure will be as important as gravitational 
perturbations. The purpose of this task was to determine the 
effects of solar radiation pressure on the SI’S orbit. It 
was shown that the eccentricity of t lie orbit can get rather 
large (.08) even when it is Initially /,ero. This is the pri- 
mary difference between the SPS orbit and other geosynclironous 
sat el 1 i te orbits. 

The SPS configuration being considered here is described 
in a study report by the Johnson Space Center. However, the 
results in this study are applicable to any geosynchronous 
satellite that resembles a flat surface that continually faces 
the sun. The main purpose of this study was to investigate 
the orbital evolution of the SPS over its expected thirty year 
lifetime. An additional goal was to describe the mf)tion with 
analytical formulas. These formulas could then be used as a 
basis for developing an orbit control theory that will mini- 
mize station keeping costs. 

The perturbing forces acting on the satellite were studied 
in detail. To a first approximation, three types of forces 
can be considered separately since they have different effects 
on the orbit . 

(1) Ix)ngitude dependent tesseral terms in the earth's 
geopotential cause a slow drift of the satellite's 
mean longitude. 
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(2) Sun and moon Kravlty cause a rotation of the 
orbital plane. 

(3) Solar radiation pressure will cause an incrc*ase 
in the orbital eccentricity. 

The following is a summary of the main results of the 
study. The dc'tailed rejK)rts are given in references 8 and 9, 

(1) An analytical solution was developed, giving the 
motion of orbital eccentricity and lineof apsides as a func- 
tion of lime. This solution has a concise formulation and 
is accurate to within a few percent. 

(2) It was shown that the eccentricity of the sun's 
orbit about the earth causes a linear increase in the eccen- 
tricity of the SPS orbit. 

(3) The eccentricity and inclination of the SPS orbit 
can be control lc*d by choosing certain nonzero values of ec- 
centricity and inclination. 

(4) Graphical data on the orbital elc‘ments was put into 
appropriate form for use in mission planning or feasibility 
studies. 

(5) Explicit equations were develop<*d for the dally 
motion in longitude and latitude. It was shown that, for 
small eccentricities and inclinations, the ground track will 
be an ellipse centered on the equator and mean longitude. 
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6.0 A HKKINKD Nt MKHICAL MKTHOD FOR CAf>CULATIONS 
OF SATKLLITt ORBITS 

Solutions to ordinary differential equations are usually 
obtained through the application of either analytical or nu- 
merical solution methods. Each approach has its own set of 
advantages and disadvantages. A particular method Is usually 
chosen according to how it may satisfy foreseen applications. 

Analytl<*al methods require extensive development of math- 
ematical formulas. Once they have been obtained, the formulas 
can be usc‘d to study certain global properties of the solu- 
tion. Also, the analytical formulas can be used in computing 
machines to provide extremely rapid numerical calculations. 

There is relatively little apriori mathematical develop- 
ment of numerical methods. They allow extremely accurate nu- 
merical calculations of the solution and are usually not dif- 
ficult to program on the computer. However, the solution must 
be developed step-by-step, which can take large amounts of 
computer time. Also, very little qualitative or global infor- 
mation on the solution is available. 

For some applications, it may be advantageous to develop 
new specialized methods that combine features of both the ana- 
lytical and the numerical methods. This is the suggested ap- 
proach for those applications that require repetitive solutions 
of the differential equations, but for which no analytical so- 
lutions are available. Therefore, the advantage in speed of 
an analytical solution might be partially realized by using 
a limited application of an analytical method. Numerical step- 
by-step calculations still need to be done, but usually in a 
more efficient manner. 

This task concerns a new approach for the solution of ar- 
tificial satellite trajectory problt'ms. The basic idea is to 
apply an analytical solution method (the Method of Averages) 
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to an appropriate formulation of the orbital mechanics equa- 
tions of motion (the KS-element differential equations). The 
result is a set of transformed equations of motion that are 
more amenable to numerical solution. 

An important problem that is found when attempting the 
numerical solution of the satellite differential equations, 
is that the maximum si/.e of the numerical integration steps 
is limited by the high frequency (short period) tesseral terms 
contained in the geopotential. This problem appears even with 
the KS formulation. Even though their amplitudes are small 
and may be considered negligible for many applications, the 
short period terms cause the following practical probl<*ms: 

(i) Since the tt*sseral terms depend on time, the orbital 
frequency bast'd on the total energy is no longer 
constant. If small steps are not taken to "track" 
the high frequency oscillations, large down range 
errors can result. 

(ii) If the tesseral terms were simply neglected in the 
numerical integration, larger steps could bo taken 
but unacceptable intrack errors would result. 

(Hi) Evaluation of all tesseral terms in the numerical 
integration force model consumes a major part of 
the computation time. 

An approach to solving the three above problems is de- 
veloped in this task. Basically, the idea is to carry out a 
numerical integration without the tesseral terms. In order 
to avoid the second of the above problems, the KS elements 
are initialized with a mean frequency (total energy). The 
mean values are obtained via a specialized application of the 
Method of Averages. Since the force model contains only zonal 
geopotential terms plus, possibly, atmosphei'ic drag and ex- 
ternal bodies, computer run time and stability problems asso- 
ciated with the tesseral terms are avoided. 

A Modi‘'<t'd Method of Averages was developed and applied 


to the KS d i f ferf>nt ial equations. This theory was then used 
to develop an Averaginfc Initialisation Algorithm. Ti*stH show 
that this algorithm can decrease the computation time required 
for prc'dicting near earth orbits. For pr<-diction intervals 
of one or more days, the computation time will be reduced by 
a factor of three to four. The mc*thod is particularly suited 
to iterative applications whore accurate in-orl)it satellite 
position information is reejuired. 

Documentation of theoretical developments and numerical 
experiments is given in reference 10. The computer subroutine 
package for the initialisation algorithm has been Implemented 
in the KSFAST program and is documented in reference 13. 
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7.0 A COMPARATIVE STUDY OF THE UNIFIED SYSTEM FOR 
ORBIT COMPUTATION AND THE FLIGHT DESIGN SYSTEM 


The high flight rate anticipated for the Shuttle Trans- 
portation System (STS) era beginning in the 1980's is caus- 
ing a review of the mission design and analysis computer soft- 
ware structures. The rather individual case-by-case planning 
procedures and software products which were us<?d for earlier 
missions are no longer ado(|uate when considering the project<*d 
budgets for the STS operations. Standardized planning tools 
are absolutely necessary in the environment of shuttle flight 
design since many tasks are rc'petitive. Also, profit must be 
taken from the experience gained with earlier missions and from 
the advances made in the areas of flight design techniques 
and computer hardware and software standards. 


In this study, two recent software systems were described 

cur- 

MPAD of NASA/JSC'*’ . The other 

cur- 

MAD 

report on the results of this study is given in reference 11 


and compared. One is the Flight Design System (FDS) 
rently being developed at MPAD of NASA/JSC^ 
is the Unified System for Orbit Computation (USOC) 
rently operational at MAD of ESA/ESOC . The complete 


The FDS software consists of a set of functional and 
utility processors, components of the data base structure and 
the executive logic. It is designed as an ef f Ic lent production 
tool for flight design and analysis tasks, to be accomplished 
for the support of operational shuttle missions and orbital 
flight tests. It also includes an interface with a system for 
automated generation of flight planning documents. 


Mission Planning and Analysis Division of NASA/Johnson 
Space Center, Houston, Texas. 

Mission Analysis Division of the European Space Agency/ 
European Space Operations Centre, Darmstadt, Germany. 
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USOC has been UHed for the planning, analysiR and doKl^n 
of Hpacecraf t mlHsionH, both earth orbital and interplanetary. 
It tnahes use of an automated approach to the ^gneration of 
probl<*m specific application programs f r«>m a library of func- 
tional modules. A mtjdule selector program interfaces with 
the user and selects the necessary modules for a particular 
application. The USOC design expressly facilitates (juick 
m<idi f ications by the analyst. 

The FDS is a produotion orienlc*d software product, 
mainly designed to satisfy flight design needs of the class 
of shuttle missions which are considered as standard missions 
or deviations of such. It is not intended to satisfy the total 
analysis nec*ds of MPAD . One of the goals of this study was 
to evaluate the approach of USOC for the planning and design 
of unique shuttle flight phases. 

The two software systems were compared on the basis of 
the following considerations: 

• Scope of applications 

• Software system structures 

• Hardware/software requirements to host computer system 

• Software flexibility 

• System transportability 

• Unique features of each system 

Although there exists some overlapping of the scope of 
applications of USOC and FDS , Uje structures of these two 
systems are very different, according to their different basic 
requirements. Their requiremt-nts are in fact, incompatible. 
While the FDS is designed to meet the requirements of a 
standardized production tool, USOC is designed for rapid 
generation of particular application programs. The main em- 
phasis of USOC is put on the adaptability to new types of 
missions . 


A software system having a USOC-like structure, adapted 
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to the specific m‘eds of MI’AD , would be appropriate to 
support planning tasks in the area of unique STS missions. 
There appears to be a need for such an additir)nal system 
within MPAD 


8.0 ANALYTICAL FORMULATION OF SELKCTKD ACTIVITIES 
OF THE REMOTE MANIPULATOR SYSTEM 


The Remote Manipulator SyHt<^m (RMS) will bo UHori by 
NASA for Shuttle on-orbit payload (P/L) iloyment and re- 
trieval. Lart;e Himulation software is beii.». used by Boveral 
orK»n l 2 at ions to support the development of the RMS . How- 
ever, for certain analyses such as evaluation of optimum RMS 
maneuvers, interactions with attitude and orbit control of 
the shuttle, etc., much simpler and faster mathematical models 
are needed. In addition, analytical closed form solutions 
will allow an estimation of operational limits of the system. 
These limits are not easily detectable with the big simula- 
tion programs that require considerable computing time. Fi- 
nally. the transparency and redundancy of a simplified math- 
ematical model are further advantages to be expected. 

The purpose of this task was to develop new meth<.ds of 
analysis for the Shuttle Orblter basc*d RMS . The* methods 
are basi-d on simplified kincmiatlcal models. The motions of 
the simplified mathematical model are described either by ex- 
pl Iclt analytical formulas or by concise numerical algorithms. 

The purposes of the new methods are to: 

• Provide insight into the coupled motion of the 
Orblter - RMS - Payload system. 

• Provide fast answers to typical problems in the 
design and analysis of Shuttle missions requiring 
the RMS. 

• Provide qualitative and quantitative results on 
the motion to be used as a check on numerical 
simulation software. 

• Isolate and analyze those RMS operations that 
are most critical from the point of view of 
mission success or RMS design limits. 
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• Invf'Htiga. ;'irthc*r fx* of th** •«»., he- 
mal leal model in onler ' * ,,vc account oi 

arm f lox Ibl 1 1 ty , drive motor characterlnt Ich 
and KMS control lawn. 

The equal lonH of motion for a mineral treatment of the 
couplt*d iTK)tlon of Orblter and P/b were derived. It was shown 
that P/L capture and P/L deployment are Koverm>d by the 
same equations. Also, all phaHos of coupled motion are de- 
scribed by the same equation of motion. 

In order to demonstrate the applicability and t»lficiency 
of the derived equations, some example problems are given. 

The corresponding equations were simplified to the extent of 
these special cases. Explicit analytical solutions are pre- 
sented for the most interesting HMS phases of operation. 

These are the following: 

• P/L capture, 

• P/L stowage into cargo bay and its removal frt>m 
the cargo bay, 

• P/L deployment. 

The analytical formulations were used to study problems 
in mission planning and on-orbit operations. Recommendations 
are given on maximum capture velocity of payloads, optimum 
capture philosophy, payload stowage and deployment techniques, 
and hardware modifications. A dc’tailed report on the study 
is given in reference 12. 






9.0 MAINTKNANCE AND MODIFICATION 
OF THE KSFAST PROGRAM 


Since the initial reloaHe of the KSFAST program, a few 
minor problems have been detected, recommendat iona for improvo- 
mentH have been made and new theories for the KS elements 
have been developed and proven. Therefore, the KSFAST pro- 
gram has been modified to Incorporate some of these imprt>vo- 
inents and to correct deficiencies. 

There have been three basic areas of modification: 

a) Changing some of the numerical constants to the 
more generally accept«*d values. 

b) Changing the input/output units from Earth radii 
and minutes to Kilometers, soconds and lays. 

c) Incorporation of the mean element thc?ory for 
approximating some of the Earth's geopotential 
terms . 

The first two modifications were done to make the KSFAST 
program more acceptable by the using public. The third in- 
volved a considerable programming effort to implement, bat 
resulted in significant savings in program execution time. 

(For a complete description of the mean element theory, see 
reference 10). 

A revised version of the KSFAST program documentation 
(reference 13) has been published, including all recent modl- 
f icat ions . 
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